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Abstract: Regional tree-ring width chronology of the Scots pine (Pinus sylvestris L.) was constructed from 
8 sites in the forest-steppe belt situated in the foothills of the Selenga River basin, Russia. Moisture 
information contained in tree-ring width chronology was obtained through linear regression 
reconstruction models of annual August—July precipitation and annual water discharge of the Selenga 
River during the period 1767—2015. Comparison of the smoothed series allowed estimating long-term 
variation component of these moisture regime parameters with a high precision. At the same time, 
regional drought indices are less correlated with pine radial growth, because they have contribution of the 
other environmental variables, which are much less reflected in the tree-ring of the investigated pine forest 
stands. Reconstructed dynamic of the moisture regime parameters is supported by documental evident of 
many socially significant events in the regional history, such as crop failures caused by both droughts and 
floods, and catastrophic fire in the Irkutsk City in 1879. Also, dependence of the amount of precipitation 
in the study area on the atmospheric circulation in Central Asia is revealed to have a similar pattern with 
other regions, i.e., a negative relationship of precipitation with the development of large high atmospheric 
pressure area within its center in the Altai and Tianshan mountains. 
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1 Introduction 


Recent substantial climatic changes support great scientific interest in the analysis of past climate 
dynamics. However, main climatic variables have distinct characteristics. Indeed, numerous 
studies have shown that temperature variation has a high spatial coherence up to the global level 
(Easterling et al., 2000; IPCC, 2014; Wu et al., 2014). In contrast, precipitation dynamics under 
the influences of relief irregularities and atmospheric masses circulation have a significant spatial 
heterogeneity, especially in mountain systems and adjacent regions. Thus, analysis of its 
dynamics is more expedient at the regional scale (Pederson et al., 2001; Batima et al., 2005; 
Dulamsuren et al., 2010). 

High latitudes and altitudes are usually characterized by sufficient humidity, and 
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spring-summer temperatures are the main limiting factor for vegetation growth. Therefore, 
contemporary warming trends lead to the positive response of ecosystems there (Wieser and 
Tausz, 2007; Bellard et al., 2012; Huhtamaa and Helama, 2017). On the other hand, there is a 
wide belt of semi-arid and arid regions in the temperate latitudes, where vegetation productivity is 
limited by the amount of available water (Sun and Liu, 2014; Yadav et al., 2015). Moisture regime 
there stems from interaction of temperature as the main evapotranspiration regulator with 
precipitation and hydrological network as sources of water (De Girolamo et al., 2011; Wang et al., 
2016). An example of such a region is the downstream basin of the Selenga River, located in the 
western Transbaikalia, Siberia, Russia. It is of great research and practical interest, because 
Selenga River provides more than half of the inflow in the Baikal Lake—one of the greatest 
freshwater lakes and a UNESCO (United Nations Educational, Scientific, and Cultural 
Organization) World Heritage Site (Berezhnykh et al., 2012; Törnqvist et al., 2014). It has led to 
the scientific attention to the hydrological and climatic regime of the Selenga River basin 
(Törnqvist et al., 2014; Frolova et al., 2017a, b; Kasimov et al., 2017). Their current change is 
evident, and its consequences are substantial not only in aquatic ecosystems (Khazheeva and 
Plyusnin, 2016; Moreido and Kalugin, 2017), but also in terrestrial ones, particularly forests 
(Mikhailova et al., 2013; Puntsukova et al., 2015). 

In consequence of the limited availability and time span of instrumental climatic and 
hydrological data series, analysis of their past dynamics often includes various natural proxy 
records (Bradley, 1999; Leclercq and Oerlemans, 2012; IPCC, 2014). Tree rings have some 
advantages, such as annual temporal resolution, precise dating, and widespread distribution of 
woody plants (Cook and Kairiukstis, 1990). Besides that, trees as the basement of many terrestrial 
ecosystems reflect the dynamics of their status and productivity (Babst et al., 2014; Klesse et al., 
2016). 

In the downstream basin of the Selenga River, dendroecological studies (Andreev, 1999, 2001a, 
b) were recently continued by our group (Demina et al., 2017). Also tree-ring-based analysis of 
climate dynamics was performed in the river upstream basin (Davi et al., 2015). Vegetation of the 
downstream basin of the Selenga River is represented by steppes and forests depending on the 
elevation. Trees growing on the border between these ecosystems (i.e., forest-steppe ecotone) are 
the most promising indicators of the moisture supply. It was of interest to analyze regional 
moisture regime dynamics using instrumental series of climatic and hydrological variables, as 
well as their reconstruction models based on tree-ring chronologies. 


2 Materials and methods 


2.1 Study area 


The study was performed in the Selenga River basin, southeast from the Baikal Lake, Russia (Fig. 
1). Its valleys are covered by steppes and surrounded by the mountains up to 1300-1800 m a.s.l., 
and the mountains are covered by Scots pine (Pinus sylvestris L.) and Siberian larch (Larix 
sibirica Ledeb.) forests. Soils are represented by mountain chernozems and mountain permafrost 
taiga turf soils. There are also aeolian sand deposits along the valleys. 

The climate is temperate continental (Alisov, 1956). It is characterized by large daily and 
annual temperature variation amplitudes with a long low-snow winter and relatively hot summer. 
January temperature ranges from —26°C to —20°C; and winter (November—March period of 
negative temperatures) precipitation is 20-35 mm. July temperature ranges from 16°C to 20°C; 
and summer (June—August period of temperatures above 15°C) precipitation ranges from 170 to 
230 mm with the maximum in July (66—90 mm). The annual mean temperature ranges from —3°C 
to 0°C, and mean annual precipitation is 240-380 mm. 


2.2 Data sources 


Monthly data series on precipitation (P) and mean temperatures (T) were obtained from Kyakhta 
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(KH, #30925; 50.4°N, 106.4°E; 791 m a.s.l.; 1940-2015) and Ulan-Ude (UU, #30823; 51.8°N, 
107.6°E; 514 m a.s.l.; 1921-2015) state weather stations in the study area; daily P and T series 
(1921-2015) were obtained from Ulan-Ude station. These stations are situated in forest-steppe 
zone and have the longest monthly series in the study area, which are sufficiently correlated 
between each other. Regional monthly P and T series during the period 1921—2015 (Fig. 1b) were 
computed from the stations data by the method of averaged standardized anomalies (Jones and 
Hulme, 1996). For comparison monthly series of P, Palmer Drought Severity Index (PDSD and 
Standardized Precipitation-Evapotranspiration Index (SPEI) during the period 1901-2015 were 
obtained from Climate Research Unit Time-Series (CRU TS 4.01) database 
(http://climexp.knmi.nl) for the area of 5O°N—52°N and 106°E-109°E. Other two drought indices 
were calculated from regional monthly P and T series from May to September, i.e., Selyaninov 
hydrothermal coefficient (HTC) and wetness index (WI) (Selyaninov, 1937; Lei et al., 2014). As a 
main hydrological characteristic of the Selenga River, its mean annual discharge (Q) from the 
Mostovoy state hydrological station (Ulan-Ude City) during the period 1934-2015 was used. 

Wood samples of pine were collected from 8 sites (Fig. 1a; Table 1) in the forest-steppe zone of 
basin, where trees are expected to be moisture-sensitive. The forest-steppe is characterized by 
low-productive open pine dry stands on sandy soils. The sparse underbrush is represented by 
Dasiphora fruticosa L., Rosa acicularis Lindl., Spiraea trilobata L., and Cotoneaster 
melanocarpus Fisch. ex Blytt. The grass has a total coverage of 5%-15% and predominantly 
consists of steppe species, i.e., Stipa krylovii Roshev., Koeleria macrantha (Ledeb.) Schult., 
Scabiosa comosa Fisch. ex Roem. & Schult., and Thymus serpyllum L. The first sample collection 
was carried out in 1999 by Andreev et al. (1999). To prolong series, we performed an additional 
sample collection in 2015, but not at all sampling sites (Demina et al., 2017). Both times, not 
suppressed healthy living trees were selected, and cores were sampled with the 5-mm-diameter 
increment borers at breast height (~1.3 m). 
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Fig. 1 Satellite map (a) with marked sampling sites (triangles), weather stations (circle) and hydrological station 
(shaded circle); location of the study area (b); climatic diagrams of 1921—2015 regional monthly precipitation (P, 
cl) and temperature (T, c2). Horizontal lines are medians, boxes are lower and upper quartiles, whiskers represent 
dataset without outliers and markers are outliers. 
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Table 1 Sampling sites and statistics of the residual tree-ring width (TRW) chronologies 


Geographical coordinates Sample Chronology EPS>0.85 period 
= Latitude Longitude oa Period ar Aba aee stdev sens r-bar Period patie 
HUH 50.9°N  106.1°E 700 1768-2015 248 55/73 0.34 0.43 0.60 1770-2014 4 
HUR 50.7°N 106.5°E 850 1751-2015 265 32/48 0.33 0.40 0.56 1777-2015 5 
KEY 51.7°N 107.2°E 700 1807-2015 209 25/42 0.35 0.36 0.52 1818-2014 6 
PDL 51.0°N 107.1°E 600 1760-2014 255 36/51 0.34 0.42 0.57 1817-2014 4 
ZGR 51.8°N 108.3°E 600 1645-2014 370 27/42 0.27 0.32 0.58 1807-2014 T 
SHN 52.1°N 108.9°E 571 1666-1998 333 11/11 0.26 0.30 0.49 1834-1998 7 
TGN 51.5°N 106.9°E 747 1811-1998 188 13/13 0.40 0.46 0.62 1837-1998 4 
ZAG 51.9°N 108.2°E 675 1719-1998 280 18/18 0.39 0.49 0.70 1726-1998 3 
Regional - - - 1645-2015 371 217/298 0.47 0.29 0.39 1769-2015 22 


Note: stdev, standard deviation; sens, mean coefficient of sensitivity; r-bar, mean inter-series correlation coefficient, EPS, expressed 
population signal. "-" means no value. 


2.3 Processing of tree-ring data 


The processing of wood samples was performed using standard dendrochronological methods 
(Cook and Kairiukstis, 1990). TRW was measured to the nearest 0.01 mm using LINTAB and 
TSAP Win program (Rinn, 2003). Cross-dating of the TRW series and measurement accuracy 
tests were performed by COFECHA program (Holmes, 1983; Grissino-Mayer, 2001). Each 
individual series was detrended using a cubic smoothing spline with a frequency response of 0.50 
at 67% of the series length as an estimation of the age trend (Cook and Kairiukstis, 1990). Then 
autocorrelation was removed from each series, because unlike TRW climatic and hydrological 
variables in the study area have not any significant autocorrelation component. Resulting residual 
individual series contain only interannual variability (Vaganov et al., 1996). On their base local 
TRW chronologies were obtained for all sites as a bi-weighted mean (Cook, 1985). To analyze 
coherence among local chronologies in different frequency domains, we decomposed these data 
into high- and low-frequency components using binomially weighted 5-years reciprocal filters 
(see examples of 5-, 8- and 13-year reciprocal filters using in LaMarche and Fritts, 1972; 
Schweingruber, 1988; Yuan et al., 2013; Zhang et al., 2014). As correlations between local 
chronologies are significant in all frequency domains (Table 2), we were able to combine them in 
one regional chronology. 

The following statistical characteristics of TRW chronologies were used, i.e., standard 
deviation (stdev) as estimation of general variability of TRW; mean coefficient of sensitivity (sens) 
to the environmental interannual variation; and mean inter-series correlation coefficient (r-bar) as 
estimation of variation consistency between trees (Cook, 1985). The expressed population signal 
(EPS) was used to identify the period when chronology contains enough samples to capture 
sufficient percentage of the external signal common for the population on local or regional scales 
(Wigley et al., 1984). 


2.4 Dendroecological analysis and reconstruction 


Monthly P and T series in the study area are significantly negatively correlated (r= —0.38— —0.17) 
during summer. It is a typical pattern for inter-mountain valleys of the South Siberia (Bazhenova 
and Tyumentseva, 2010). Therefore, to estimate climatic response of the tree growth, we 
calculated paired and partial correlations of TRW chronology with P and T respectively with 
Seascorr program (Meko et al., 2011). It also integrated climatic series for seasons of various 
lengths. Seasons of 1, 3, 6, and 12 months were chosen in this study. 

In correlation analysis of environmental and tree growth variables, cross-correlation with time 
shifts (lags) was used to take into account possible temporal differences between variables. To 
obtain tree-ring-based reconstructions the moisture regime parameters, we used linear regression 
models. The model fitting to instrumental data was estimated with coefficient of determination 
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Table 2 Correlation coefficients between unfiltered, high-, and low-pass filtered TRW chronologies during the 
period 1813—1996 


ZGR SHN ZAG HUR PDL HUH KEY TGN 
Unfiltered 
SHN 0.87 
ZAG 0.71 0.68 
HUR 0.49 0.46 0.49 
PDL 0.54 0.50 0.51 0.76 
HUH 0.45 0.41 0.45 0.79 0.66 
KEY 0.58 0.56 0.55 0.54 0.57 0.61 
TGN 0.39 0.35 0.53 0.63 0.56 0.63 0.64 
Regional 0.67 0.63 0.69 0.88 0.81 0.88 0.73 0.74 
Low-pass filtered 
SHN 0.88 
ZAG 0.78 0.70 
HUR 0.47 0.42 0.47 
PDL 0.60 0.53 0.63 0.76 
HUH 0.47 0.34 0.48 0.78 0.73 
KEY 0.65 0.52 0.55 0.55 0.60 0.61 
TGN 0.38 0.25 0.47 0.60 0.59 0.61 0.54 
Regional 0.69 0.59 0.70 0.86 0.85 0.86 0.71 0.71 
High-pass filtered 
SHN 0.89 
ZAG 0.70 0.70 
HUR 0.53 0.52 0.49 
PDL 0.51 0.50 0.42 0.75 
HUH 0.45 0.49 0.47 0.78 0.60 
KEY 0.57 0.62 0.58 0.54 0.54 0.61 
TGN 0.48 0.49 0.60 0.64 0.54 0.64 0.70 
Regional 0.69 0.70 0.70 0.89 0.78 0.88 0.75 0.78 


Note: All correlations are significant at P<0.05 level. 


(R?), adjusted coefficient of determination (R?aaj) and F-test. The stability of model was analyzed 
using split-sample procedure (Meko and Graibill, 1995). The total length of instrumental series 
was divided into two sub-periods. The model was calibrated on one of sub-periods with the value 
of R?aj, validation was estimated with correlation coefficient (r) between model and instrumental 
data for the other sub-period, reduction of error (RE), and coefficient of efficiency (CE) statistics 
(Cook et al., 1994). To analyze the frequency domains of variation for instrumental and 
reconstructed series, we deconstructed with aforementioned 5-year reciprocal filters. Then we 
performed sign test (ST) for full series and their high- and low-frequency components. It is a ratio 
of agreements to disagreements between instrumental and reconstructed series in signs of 
departures from their mean values over considered period (Cook and Kairiukstis, 1990; Buckley, 
2009). 

The influence of the atmosphere circulation on the regional moisture regime and on TRW was 
investigated using composite map of the 500-mb geopotential height anomalies. It was calculated 
as differences between extremely dry and extremely wet years during the period 1948-2015 
(period of the reanalysis database NCEP/NCAR) with online mapping tool from the NOAA Earth 
System Research Laboratory (https://www.esrl.noaa.gov/psd/cgi-bin/data/composites/ 


printpage.pl). 
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3 Results 


3.1 TRW chronology development 


The length of local Pinus sylvestris TRW chronologies ranges from 188 (TGN) to 370 years 
(ZGR). Statistical characteristics of local chronologies (Table 1) are acceptable for the 
dendroecological research. All chronologies are significantly and positively correlated (Table 2) 
and well cross-dated, which justifies possibility of their generalization into one regional 
chronology. Regional chronology has lower values of sensitivity (sens=0.29) and inter-series 
correlation (r-bar=0.39) than local ones (sens=0.32—0.49; r-bar=0.49-0.70). On the other hand, 
standard deviation value is higher for regional chronology than for local ones (stdev=0.47 and 
stdev=0.26—0.40, respectively). However, all regional chronology characteristics are acceptable as 
well. 


3.2 Climatic response of TRW and reconstruction of annual precipitation 


Seascorr analysis showed that regional TRW chronology has a positive response to P. It is 
maximal in June of current year and August-September (Fig. 2a) of preceding year. Partial 
correlations showed that T influences pine growth independently from P only in winter 
(October—February), and this impact is positive. Generalization for longer seasons revealed 
maximal correlations of TRW with P for July—June and August—July annual periods. To clarify 
reason behind approximate equality of these periods and to choose one of them for the following 
analysis, we calculated annual P series from daily data for periods ending 30 June—31 July with 
1-day intervals of moving. Correlations of these series with TRW have maximum for periods 
ending 18—21 July (Fig. 2b), so actual influence of precipitation on the development of single 
tree-ring of pine in the study area does not coincide with the calendar months but lasts 
approximately from the third decade of preceding July to the second decade of current July. 
Therefore, P reconstruction was performed for August—July (Fig. 3). 

Linear regression reconstruction function is P=142.5+164.3xTRW (Figs. 3al and a2). Its 
statistical characteristics and results of split-sample validation procedure are shown in Table 3. 
Sign test showed that reconstruction model captures precipitation variation in both high- and 
low-frequency domains. 

Grid series of August—July P for 50°N-52°N, 106°E-109°E area have correlations with 
reconstructed series (r=0.59—0.69), and instrumental series (r=0.70-—0.85) . 
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Fig. 2 Climatic signal in regional tree-ring width (TRW) chronology. Seascorr summary (a): paired correlations 
(rxp) between regional chronology and precipitation for periods of 1, 3, 6, and 12 months lengths (upper panel); 
partial correlations (rxr) between regional chronology and temperature for periods of 1, 3, and 6 months lengths 
(lower panel). Asterisks mark the months of preceding year. Correlation between regional chronology and P 
summarized for moving annual periods from 1 July—30 June to 1 August-31 July with 1-day intervals (b). Periods 
with maximal correlations are marked. 
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Fig. 3 Dynamics of the regional moisture regime. Regional August-July precipitation (P) instrumental series 
and TRW-based reconstruction model (al, a2). Low-frequency components of the mean annual discharge (Q), 
August—July P (b) and regional tree-ring width chronology (TRW, c) of the Selenga River, estimated with 5-year 
bi-weighted low-pass filter. On the panel (al and a2) horizontal solid line is the mean and dashed lines are the 25* 
and 75" percentiles of the regional instrumental P during the period 1922-2015. Uncertainty in reconstructed 
values is shown by an 80% confidence interval (shade). 


Table 3 Statistical parameters of split-sample calibration—verification for tree-ring-based reconstruction of 
August—July precipitation 


Calibration Verification 


ST ST 
Period R Rg F Period r RE CE 
Full High Low Full High Low 


1922-1968 0.30 0.28 19.1 36/11 28/17 38/7 1922-1968 0.55 0.35 0.27 34/13 35/10 38/7 
1969-2015 0.63 0.62 75.5 34/13 26/19 38/7 1969-2015 0.79 0.64 0.60 34/13 34/11 38/7 
1922-2015 0.44 0.44 72.9 71/23 65/25 76/14 


Note: R?, coefficient of determination (explained variance); Roads adjusted coefficient of determination; F, Fisher test statistics; ST, sign 


test, calculated for unfiltered (full), high-, and low-pass filtered series; r, Pearson correlation coefficient; RE, reduction of error; CE, 
coefficient of efficiency. “All F-test values are significant at P<0.05 level. 


3.3 Historical documentation of extreme years 


Many extremes of both signs of the reconstructed P outside the instrumental data period are 
confirmed by historical sources (Table 4). For example, some maxima are consistent with the 
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reports on floods from high water or heavy rains in Transbaikalia and Irkutsk provinces. 
Kozhuhov (1967) mentioned that "1801, 1811, 1820, 1830, 1840, and 1848 were the most severe 
years for agriculture". There is also stated that in 1802, in Transbaikalia "due to lack of bread in 
the summer time, people were fed with the quinoa and other grasses"; in 1809 "the settlers 
suffered a calamity of bad harvest" which lasted several years; "crop failure continued for several 
years and almost all Transbaikalia suffered from it in 1811 and during the period 1813-1816"; 
"from the bad harvest of grass and the lack of a forage feed in September 1838, cattle began to 
die...". Poor harvests of 1811—1816 were also mentioned by Shcheglov (1884). As reconstructed P 
showed, bad harvest could have been caused by drought, as well as by heavy spring-summer rains. 
Besides that, catastrophic fire of 1879 in Irkutsk was also drought-related (Kudryavtsev and 
Vendrih, 1971). May and June of that year were extremely hot and dry, and when fire ignited (22 
and 24 June) the weather was also very hot, which made fire-fighting difficult. As a result, the 
fires produced a huge devastation, more than half of the city was destroyed (Antonov, 1996). 


Table 4 Comparison of reconstructed precipitation (P) extremes with data on historical events 


Year P extremes Historical event Reference 
1785 + Heavy rains, flood Pezhemsky and Krotov (1911) 
1801 — Severe year for agriculture Kozhuhov (1967) 
1802 - Food deficit Kozhuhov (1967) 
1806 + Heavy rains, flood Pezhemsky and Krotov (1911) 
1809 - Bad harvest Kozhuhov (1967) 
1810 + Bad harvest Kozhuhov (1967) 
1811 — Severe year for agriculture Kozhuhov (1967); Shcheglov (1884) 
1812 + Bad harvest Kozhuhov (1967); Shcheglov (1884) 
1813 + Bad harvest Kozhuhov (1967); Shcheglov (1884) 
1814 + Bad harvest Kozhuhov (1967); Shcheglov (1884) 
1815 — Bad harvest Kozhuhov (1967); Shcheglov (1884) 
isie P Bad harvest Kozhuhov (1967); Shcheglov (1884) 
Heavy rains, flood Kozhuhov (1967) 
iss m Heavy rains, flood Shcheglov (1884) 
Severe year for agriculture Kozhuhov (1967) 
1827 + Heavy rains, flood Pezhemsky and Krotov (1911) 
1828 + Heavy rains, flood Pezhemsky and Krotov (1911) 
1829 + Heavy rains, flood Pezhemsky and Krotov (1911) 
iss i Heavy rains, flood Pezhemsky and Krotov (1911) 
Severe year for agriculture Kozhuhov (1967) 
1838 - Bad harvest Kozhuhov (1967) 
1840 7 Severe year for agriculture Kozhuhov (1967) 
1841 + Flood Kozhuhov (1967) 
1848 — Severe year for agriculture Kozhuhov (1967) 
1879 — Catastrophic fire in Irkutsk Kudryavtsev and Vendrih (1971); Antonov (1996) 
1897 + Heavy rains, flood Voznesensky and Shostakovich (1913) 


Note: Plus (+) represents positive P extremes, and minus (—) represents negative ones. 


3.4 Anomaly of atmospheric circulations in extreme years 


We chose extremely dry 1969, 1980, and 1987 and extremely wet 1962, 1985, and 1991 to 
analyze the effect of atmospheric circulation on the regional moisture regime. These years have 
extremal amount of precipitation in April—July, as well as extremes of the same sign in the 
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regional TRW chronology (values outside the range of mean+stdev). As shown on the composite 
500-mb geopotential height anomaly map for the same months (Fig. 4), drought is associated with 
the high atmospheric pressure area with center in the Altai and Tianshan mountains. 
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Fig. 4 Composite map of April-July 500-mb geopotential height difference between dry (1969, 1980 and 1987) 
and wet (1962, 1985 and 1991) years. Color bar represents departures (m) from 1981-2010 climatology. Circle 
indicates study area location. 


3.5 Relationships of precipitation and its reconstruction model with drought indices 


Four indices describing the moisture availability on the territory were considered. HTC and WI 
can be calculated only for a period of temperature above +10°C, therefore in the study area they 
were obtained from May to September. HTC is persistently more closely related to August-July P 
than WI. The correlations of these indices with instrumental P series are r=0.61—0.78 in current 
May-July and 7=0.54—0.65 in preceding August-September, and the best convergence with the 
annual precipitation is shown by the indices combined for both these periods (r=0.79—0.96). 
Comparison of these drought indices with the regional TRW chronology showed that relationships 
have the same temporal pattern during season, but they are weaker, and difference between 
indices is much smaller, i.e., r=0.43-0.45 in May-July, 7=0.40-0.47 in preceding 
August-September, and r=0.57—0.61 in the combined period. 

Spatially distributed PDSI and SPEI can be calculated throughout all the year, so for this study 
they were generalized for August—July. Nevertheless, despite more temporal accordance 
correlations of their averaged series for 5|0°N—52°N, 106°E-109°E area with instrumental P series 
are noticeably lower than correlations of the indices above, i.e., r=0.41 for PDSI and 7=0.40 for 
SPEI, the TRW chronology correlate with these indices evens less, i.e., r=0.35 for PDSI and 
r=0.27 for SPEI. However, it should be noted that all listed relationships are statistically 
significant at P<0.05 level. 


3.6 Relationships of precipitation and its reconstruction model with discharge of the 
Selenga River 


The only available series of the mean Q of the Selenga River was for calendar years, i.e., from 
January to December. Therefore, it makes sense to consider Q relationship not only with annual P 
from preceding August to current July, but also with P from current August to succeeding July 
(i.e., the same series with lag of +1-year), and to P of other annual periods in-between (Fig. 5). 
The correlations between Q and August—July P are 7=0.48 without lag and r=0.29 with +1-year 
lag, both of them are significant at P<0.05 level. The maximal correlation between Q and annual 
P is observed for October-September (7=0.65). The similar pattern is observed between Q and 
TRW, i.e., 7=0.44 without lag and r=0.32 with +1-year lag. Thus, we can describe it quantitatively 
as two-factor linear regression with the following equation: Q(TRW)=229+382xTRW 
+292xTRW.1. This regression model has R?,q=0.28 and correlation with instrumental series 
r=0.54. 
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Comparing low-frequent components of P, Q, and TRW series (Fig. 3b), they have sufficiently 
higher correlations between each other, i.e., r=0.59 between precipitation and discharge smoothed 
series, 7=0.61 between discharge and TRW ones, and r=0.79 between precipitation and TRW ones. 
Smoothed reconstructed and instrumental discharge series have a positive correlation (7=0.63) 
between each other. Also, comparison of smoothed series shows that they all have the same 
in-phase quasicyclic oscillation with period length of ~25—35 years (Balybina, 2006; Demina et 
al., 2017). 


Calendar year 
Q 


Jan 
max correlation with Q 


Oct* 


P—max dorrelation with TRW 


EBEEEAREEREREES 


ana] 


Q(TRW) linear regression 


~21 Jul* ~20 Jul 


Time 


Fig. 5 Relationships between mean annual Q (discharge) of the Selenga River, regional TRW chronology, and 
regional precipitation (P) calculated for annual periods of maximal correlations with Q and TRW. Asterisks mark 
the preceding year. For TRW chronology approximate period of maximal precipitation impact is shown according 
to daily climatic data (see Fig. 2), approximate period of TRW development is marked with dark shade (from 
T>+10° to the end of maximal P impact period). Dashed lines represent P and TRW with +1-year lag. Arrows 
represent significant at P<0.05 positive correlations. Correlation values r>0.50 are marked with bold arrows. 


4 Discussion 


The development of the regional TRW chronology made it possible to integrate the variation of 
pine growth from the entire area from which the material was collected, since the sites were 
selected in relatively close natural conditions and distributed throughout the study area (Fig. 1a). 
However, some differences between sites both in regards to habitat conditions and to local climate 
(because of uneven precipitation) lead to a decrease in the statistical characteristics of regional 
chronology relative to local ones. As a consequence, the period of its suitability for reconstruction 
also does not reach the maximum observed for ZAG local chronology. Despite that, in the 
investigation of the regional climate and the hydrological regime of territory, it is more expedient 
to use regional chronology and not local one. 

The positive reaction of chronologies on precipitation as a source of moisture is typical for the 
forest-steppe zone of South Siberia (as shown in Andreev et al., 1999; Demina et al., 2017). 
However, unlike other territories (e.g., Shah et al., 2015, Babushkina et al., 2017), the temperature 
of the warm period in the study area only indirectly affects the growth, through the regulation of 
water loss. This is indicated by the absence of significant partial correlations, while paired 
temperature—growth correlations in the study area are significant (Demina et al., 2017). A positive 
response of growth to the temperature of the late autumn and most part of the winter is due to the 
fact that in conditions of low snow cover thickness soil freezing in extremely cold winters can 
cause damage to the most active part of the fine roots of trees (Tierney et al., 2001; Schaberg et al., 
2008). Such damage can occur since October when temperature sharply decreases before 
establishment of stable snow cover, and continue until March, when severe frosts are unlikely to 
occur. 
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Clarification of the most significant period influencing by climate on tree growth using daily 
climatic series and the moving periods of their generalization was carried out in previous research 
(Belokopytova et al., 2018). The results demonstrated that climate can vary significantly between 
remote sites within the region due to climatic gradients even when their habitat conditions are 
homogeneous. Elevation and landscape features also can shift timing of cambial phenology and 
consequently tree ring development (Prislan et al., 2013; Jiang et al., 2015; Kraus et al., 2016). 
Therefore, taking into account the horizontal distance between the sites up to 170 km and the 
vertical difference up to 280 m, TRW series integrated into the regional chronology react to the 
climate during different periods even within one year. On the other hand, climatic variations lead 
to a change in the beginning and end of the tree ring development (Prislan et al., 2013; Gricar et 
al., 2014; Swidrak et al., 2014; Yang et al., 2017). Thus, the end of the maximum annual period of 
precipitation variation (18—21 July) revealed in this study gives us the estimation of the mean or 
most probable moment of the cambial activity ending for pine trees growing in the forest-steppes 
of the study area. 

For the selected August—July period, the precipitation contribution in the tree-ring development 
of pine in the forest-steppe zone is significant and sufficient to use the TRW chronology for 
precipitation reconstruction. The quality and accuracy of the reconstruction model is confirmed 
not only by its statistical characteristics, but also by convergence of reconstructed precipitation 
dynamics with other data sources, i.e., spatially distributed precipitation series and historical 
evidence of extremal events that are significant for the society and national economy of the 
region. 

Analysis of the atmospheric masses circulation revealed pattern in the spring-summer drought 
development, typical not only for the study area. Despite the choice of different dry and wet years, 
the composite map of the geopotential height has a high similarity with maps constructed for 
Khakassia (Kostyakova et al., 2018) and Kashmir, India (Shah et al., 2018). Moreover, negative 
correlations between precipitation amount and geopotential height were reported in Mongolia and 
northern China (Hu, 1997; Wang et al., 2017). Consequently, the formation of anticyclones in the 
Altai and Tianshan mountains in the spring-summer season reduces the amount of precipitation 
not only on the spot, but also spreads its influence to other regions of Central Asia from northern 
India to South Siberia, due to the flow of dry atmospheric masses there. 

It is understandable that more closely related to the regional precipitation series amongst 
drought indices are ones calculated on the basis of the same series, i.e., HTC and WI. When 
calculating WI, the precipitation variable is used as an argument of the logarithm function (Lei et 
al., 2014). This non-linearity of the relationship reduces the correlations in comparison with the 
HTC, which is proportional to the sum of precipitation. The less pronounced difference between 
HTC and WI correlations with TRW chronology indicates that precipitation during extreme 
droughts has stronger influence on the pine growth, i.e., this influence is nonlinear. The possible 
reason of this nonlinearity is that drought not only depresses radial increment rate, but also 
reduces duration of cambial activity (Rossi et al., 2014). Weaker relationships of PDSI and SPEI 
with precipitation are partially explained by their calculation on the basis of another data set, and 
partially by a greater contribution of other climatic variables to these indices, including 
temperature (Vicente-Serrano et al., 2010; Osborn et al., 2017). And since the temperature does 
not exert an independent influence on the tree growth, all four considered drought indices 
correlating with regional TRW chronology (and P reconstruction based on it) are weaker to the 
certain degree than those correlating with instrumental precipitation series. 

Similarly to the drought indices, the water balance of hydrological objects depends not only on 
the amount of precipitation, but also on evaporation, i.e., it is connected with temperature (e.g., 
water balance equations in Shanahan et al., 2007; Giadrossich et al., 2015). Another factor 
lessening the correlations of the Selenga River with regional precipitation series is that the study 
region represents only part of the Selenga basin. Therefore, discharge measured at the Mostovoy 
Station integrated precipitation amount from much larger area. As shown in Figure 5, the inflow 
of water from precipitation into the riverbed occurs with a certain lag, which can be attributed to 
two main reasons: (1) low speed of surface and underground runoff (Gillies et al., 2011, 2015; 
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Babushkina et al., 2017), and (2) solid form of winter precipitation falling which leads to delay of 
the water flow until the spring snowmelt. In connection with this lag, the average annual 
discharge of the Selenga River has significant correlations with both the August—July 
precipitation and the TRW chronology with the 1-year lag. This allowed us to develop much more 
accurate two-factor regression model of discharge reconstruction based on the tree rings 
variability with zero and 1-year lag than simple regression would be. 

A higher similarity between the low-frequency dynamics components of the precipitation, 
Selenga discharge and tree growth indicates the possibility of using the obtained reconstruction 
models to assess decadal to centennial fluctuations in the moisture regime of the territory with 
high accuracy. Such assessments are important for forecasting economic activity within the 
Selenga River basin, as well as for assessing the risks of long-term extremal variations of the 
Baikal Lake under conditions of its regulation (Törnqvist et al., 2014). 


5 Conclusions 


Radial growth of Pinus sylvestris in the forest-steppe zone of the downstream basin of the 
Selenga River is closely related to the amount of precipitation during the year, while the 
independent effect of temperature is expressed only in the winter months. Based on the regional 
TRW chronology, we reconstructed the August—July precipitation during the period 1769-2015. 
The quality and accuracy of the reconstruction model are statistically confirmed, by other sources 
of instrumental data and by documentary evidence. The historical significance for the region is 
shown not only for droughts, but also for extremely high precipitation. TRW response to other 
indicators of moisture regime, namely drought indices and Selenga discharge, is less pronounced 
than that to precipitation. The relationships of the mean annual discharge with precipitation and 
TRW were quantitatively described, resulting in a model of its reconstruction, but it has less 
accuracy than precipitation one. However, tree growth registers long-term fluctuations of both 
parameters of the moisture regime with very high accuracy. The patterns of atmospheric 
circulation in Altai and Tianshan mountains were revealed to impact the development of droughts 
in South Siberia among other regions of Central Asia. 
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